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The energy structure and the thermodynamics of ferrimagnetic Heisenberg chains of ahernating 
spins 5* and s are described in terms of the Schwinger bosons and modified spin waves. In the 
Schwinger representation, we average the local constraints on the bosons and diagonalize the Hamil- 
tonian at the Hartree-Fock level. In the Holstein-Primakoff representation, we optimize the free 
energy in two different ways introducing an additional constraint on the staggered magnetization. 
A new modified spin-wave scheme, which employs a Lagrange multiplier keeping the native energy 
structure free from temperature and thus differs from the original Takahashi Scheme, is particularly 
stressed as an excellent language to interpret one-dimensional quantum ferrimagnetism. Other types 
of one-dimensional ferrimagnets and the antiferromagnetic limit S = s are also mentioned. 



PACS numbers: 75.10.Jm, 75.50.Gg, 75.40. Cx 
I. INTRODUCTION 

Significant efforts have been devoted to syn- 
thesizing low-dimensional ferrimagnets and under- 
standing their quantum behavior in recent years. 
The first example of one-dimensional ferrimagnets, 
MnCu(S2C202)2(H20)3-4.5H20, was synthesized by 
Gleizes and Verdaguer [1] and followed by a series of or- 
dered bimetallic chain compounds [2] in an attempt to 
design molecule-based ferromagnets [3]. Caneschi et al. 
[4] demonstrated another approach to alternating-spin 
chains hybridizing manganese complexes and nitronyl 
nitroxide radicals. The inorganic-organic hybrid strat- 
egy realized more complicated alignments of mixed spins 
[5]. There also exists an attempt at stacking novel tri- 
radicals into a purely organic ferrimagnet [6]. Mono- 
spin chains can be ferrimagnetic with polymerized ex- 
change interactions. An example of such ferrimagnets 
is the ferromagnetic-antiferromagnetic bond-alternating 
copper tetramer chain compound Cu(C5H4NCl)2(N3)2 
[7]. The trimeric intertwining double-chain compound 
Ca3Cu3(P04)4 [8] is another solution to homometallic 
one-dimensional ferrimagnets, where the noncompensa- 
tion of sublattice magnetizations is of topological ori- 
gin. Besides one-dimensional ferrimagnets, metal-ion 
magnetic clusters such as [Mni20i2(CH3COO)i6(H20)4] 
[9] and [Fe8(N3C6Hi5)602(OH)i2]^+ [10], for which res- 
onant magnetization tunneling [11-14] was observed, are 
also worth mentioning as zero-dimensional ferrimagnets. 

The discovery of ordered bimetallic chain compounds 
stimulated extensive theoretical interest in (quasi-)one- 
dimcnsional quantum ferrimagnets. Early efforts [15] 
were devoted to numerically diagonalizing alternating- 
spin Heisenberg chains. Numerical diagonalization, 
combined with the Lanczos algorithm [16,17] and a 
scaling technique [18], further contributed to study- 
ing modern topics such as phase transitions of the 
Kosterlitz-Thouless type [17,19] and quantized magne- 
tization plateaux [20,21]. Alternating-spin chains were 
further investigated by density-matrix renormalization- 



group [22,23] and quantum Monte Carlo [24,25] methods 
in an attempt to illuminate dual features of ferrimag- 
netic excitations. More general mixed-spin chains were 
analyzed via the nonlinear a model [26] with particu- 
lar emphasis on the competition between massive and 
massless phases. Quasi-one-dimensional mixed-spin sys- 
tems [27,28] were also investigated in order to explain the 
inelastic-neutron-scattering findings [29,30] for the rare- 
earth nickelates i?2BaNi05. 

In order to complement numerical tools and to achieve 
further understanding of the magnetic double structure 
of ferrimagnetism, several authors have recently begun 
to construct bosonic theories of low-dimensional quan- 
tum ferrimagnets. The conventional spin-wave descrip- 
tion of the ground-state properties [22,31-33], a modi- 
fied spin- wave scheme for the low-temperature properties 
[34] , and the Schwinger-boson representation of the low- 
energy structure [35] and the thermodynamics [36] , they 
all reveal the potential of bosonic languages for various 
ferrimagnetic systems. However, considering the global 
argument and total understanding over the bosonic the- 
ory of ferromagnets and antiferromagnets [37-47], ferri- 
magnets are still undeveloped in this context especially 
in one dimension. In such circumstances, we represent 
one-dimensional Heisenberg ferrimagnets in terms of the 
Schwinger bosons and the Holstein-Primakoff spin waves. 
Based on a mean-field ansatz, the local constraints on the 
Schwinger bosons are relaxed and imposed only on the 
average. The conventional antiferromagnetic spin-wave 
formalism [48,49] is modified, on the one hand following 
the Takahashi scheme [37,38] which was originally pro- 
posed for ferromagnets, while on the other hand introduc- 
ing a slightly different new strategy [50] . The Schwinger 
bosons and the modified spin waves both interpret the 
low-energy properties fairly well identifying the ferrimag- 
netic long-range order with a Bose condensation, while 
the two languages are qualitatively distinguished in de- 
scribing the thermodynamics. We demonstrate that the 
new modified spin-wave scheme of all others depict one- 
dimensional ferrimagnetic features very well. 
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II. FORMALISM 

A practical model for one-dimensional ferrimagnets is 

two kinds of spins, S and ,s (S > s), alternating on a 
ring with antiferromagnetic exchange coupling between 
nearest neighbors, as described by the Hamiltonian, 



N 

Ti. = J ^ ] {Sn ■ S„—i + S„ ■ Sn) , 

n=l 



(2.1) 



where N is the number of unit cells. The simplest case, 
{S, s) = (1, has so far been discussed fairly well using 
the matrix-product formalism [51], a modified spin- wave 
scheme [50], the Schwinger-boson representation [36], and 
modern numerical techniques [22-25]. We make further 
explorations into higher-spin systems and develop the an- 
alytic argument in more detail. 



A. Schwinger-boson mean-field theory 

Let us describe each spin variable in terms of two kinds 
of bosons as 

4=&iAi, < = WnM-^iM), (2-2) 

where the constraints 

^ a^ana = 25 , ^ = 2s , (2.3) 

are imposed on the bosons. Then the Hamiltonian can 
be written as 

N 

n = 2NJSs -2JJ2 {^l+^n+ + ^i-^n-) , (2.4) 

n=l 

where rj„+ = (a„|6„4 - a„46„|)/2 and = 
(«nT^n-ii — ari-ii^nT)/2- The Hartree-Fock treatment 
assumes the thermal average of the short-range antifer- 
romagnetic order to be uniform and static as 

{ni^)T = {fin±)T = ^l. (2.5) 

The constraints (2.3) are correspondingly relaxed as 

N N 

^ ^ aLa„. = 2NS, ^ ^ 6^^. = 2Ns. (2.6) 

n=lo-=T,J. n=lCT=T,i 

In the momentum space the mean-field Hamiltonian 
reads 

Hmf = 2NJSs + 4Njn'^ - 4NJ{XS + us) 

— 2Jf2 ^^cosaA; {ah-^bki — (ifej.&/cT -l-H.c.) 

k 



where k is defined as nir/Na (n = 0, 1, • • • , TV — 1) with 
a being the distance between neighboring spins, and A 
and fx are the Lagrange multipliers due to the constraints 
(2.6). Via the Bogoliubov transformation 

flfcT = ak-^cosMk + 0l^smh.ek , 

<^ki = akicosh0k - f3l^smh9k , 

bk^ = /3fetCOsh6'fc - a^^sinh^fe , 

hi = /3feicosh6';o -t- a[^sinh6'fe , 



with 



tanh2^fc = 



2Q cos ak 



(2.9) 



the Hamiltonian (2.7) is diagonalized as 

Hmf = 2NJSs + ANjn^ 

-2NJX{2S -hi) - 2NJn{2s + l) + 2J^0Jk 

k 

+ -^E E {^kAa^ka+^tAjika). (2.10) 



where 



'^ka = ^k ^ wfc± (m- A); 



Wfc = V(A+7I)2^^4?)2W2^. (2.11) 
A, /i, and are determined through a set of equations 

{nl^cosh^ek + nl^smh^ek + sinh^Ok) = NS , (2.12) 

k 

J2 {nk^sinh^ek + n+^cosh^eife + sinh26'fc) = Ns , (2.13) 

fe 

E i^ka + ^ka + 1) cosh(9fesinh6lfe = NQ , (2.14) 



where the thermal distribution functions 

{al^akcr)T and n^^ = {(3l^l3ka)T are required to mini- 
mize the free energy and given by 



1 



ni.„ = 



The magnetic susceptibility is expressed as 



(2.15) 



where we have set the g-factors of spins S and s both 
equal to g. The internal energy should be given by 

E=^{Emf + 2NJSs)-2NJSs, (2.17) 

where 
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^MF = 2NJSs + - 2NJ{2XS + 2/is + A + /x) 

+2J^oJk + jJ2Yl E ^fe-'^fe^^- (2-18) 

k k r=±<j=|a 

Arovas and Auerbach [44] pointed out that relaxing the 
original constraints (2.3) into Eq. (2.6) leads to dou- 
ble counting the number of independent boson degrees of 
freedom. Therefore, in Eq. (2.17), we have corrected the 
mean-field artifact reducing the overestimated quantum 
fluctuation. 



(fc) = {S + s)cosh20k — 2\/!SscosaA;sinh20fc 

±{S~s) = LUk±iS-s), (2.26a) 

U IT 
s VS. 



± 



)- 



cosh2^fe 

cos ak sinh26fc 

(2.26b) 



B. Modified spin- wave theory: Takahashi scheme 

Next we consider a single-component bosonic represen- 
tation of each spin variable at the cost of the rotational 
symmetry. We start from the Holstein-Primakofi^ trans- 
formation 



2S 



(2.19) 



4 =bl^2s-bibn, s^„ = -s + blbn. 

Treating S and s as 0{S) = 0{s), we can expand the 
Hamiltonian with respect to 1/5 as 

n = -2NjSs + Ei + Eo + ni + no + o{s-'^), (2.20) 

where Ei and Hi give the 0{S^) quantum corrections 
to the ground-state energy and the dispersion relations, 
respectively. Via the Bogoliubov transformation 



ak = afcCosh6'fc - /5^sinh0fe 
bk = /3feCosh6'fc - Q[sinh6'fe 

they are written as 
El = -2NJ \2VSs r - (5* + s)A 

Eo = -2NJ + A^ - ( 



irA 



(2.21) 

(2.22a) 
(2.22b) 



+7i(fc)(afe/3fc + 44) 



where 



r = ^2 cos ak sinh26fc , 

k 



(2.23) 

(2.24) 
(2.25) 



71 (fc) = 2VSs cos ak cosh26'fe - {S + s)sinh26'fe , (2.27a) 

s 



70 (fc) = 



2r 



V V s 
V V s 



1'^ 



5lr-2A 



COS ak cosh26k 
sinh26»fc . (2.27b) 



The conventional spin-wave scheme naively diagonal- 
ize the Hamiltonian (2.20) and ends up with the num- 
ber of bosons diverging with increasing temperature. In 
order to suppress this thermal divergence, Takahashi 
[38] considered optimizing the bosonic distribution func- 
tions under zero magnetization and obtained an excel- 
lent description of the low-temperature thermodynam- 
ics for low-dimensional Heisenberg ferromagnets. For 
ferrimagnets, this idea is still useful [34,50] but never 
applies away from the low-temperature region as it is. 
The zero-magnetization constraint plays a role of keep- 
ing the number of bosons finite under ferromagnetic in- 
teractions but does not work so under antiferromagnetic 
interactions. Takahashi [38] and Hirsch et al. [39] pro- 
posed constraining the staggered magnetization, instead 
of the uniform magnetization, to be zero as the anti- 
ferromagnetic version of the modified spin-wave theory. 
Their scheme was applied to extensive antiferromagnets 
in both two [38,39,52-55] and one [40,41] dimensions. 
The conventional spin-wave procedure assumes that spins 
on one sublattice point predominantly up, while those on 
the other predominantly down. The modified spin-wave 
treatment restores the sublattice symmetry. We consider 
the naivest generalization of the antiferromagnetic mod- 
ified spin-wave scheme to ferrimagnets. 

The constraint of zero staggered magnetization reads 



Y^{aiar. + bibn)=N{S + s). 



(2.28) 



In order to enforce this condition, we first introduce a 
Lagrange multiplier and diagonalize the effective Hamil- 
tonian 



W = W + 2J//^(a|,a„ + 6t6„). 



(2.29) 



Then the ground-state energy and the dispersion rela- 
tions are obtained as 
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Eg = -2NJS.S + Ei; Ei=Ei+ ANJAv , (2.30) 
a;± =Si'=(A;) ; uj f (k) = ut{k) +2ucosh.2ek , (2.31) 

keeping only the bilinear terms and as 



Eg = -2NJSs + E1+E0, 



(2.32) 
(2.33) 



considering the 0(5°) interactions as well. In terms of 
the spin- wave distribution functions 



/fcBT _ I 



(2.34) 



the internal energy and the magnetic susceptibihty are 
expressed as [56] 



k T=± 



SknT 



(2.35) 



(2.36) 



k r=± 



9k, defining the BogoHubov transformation (2.21), is de- 
termined through 

71 {k) - 2r/sinh26lfc = 71 (fc) = , (2.37) 

provided we treat Ho as a perturbation to Hi. 

C. Modified spin-wave theory: A new scheme 

Although the Takahashi scheme overcomes the diffi- 
culty of sublattice magnetizations diverging thermally, 
the obtained thermodynamics is still far from satisfactory 
(see Fig. 2 later on). Within the conventional spin- wave 
theory, the quantum spin reduction, that is, the quantum 
fluctuation of the ground-state sublattice magnetization 
per unit cell, reads 



{alan)T=o = {bibn)T=o = S 
S + s 



dk 



^{S - s)"^ + 4Sssm^{k/2) 



2tt 



(2.38) 



and diverges at 5 = s. The Takahashi scheme settles 
this quantum divergence as well as the thermal diver- 
gence. However, the number of bosons does not diverge 
in the ferrimagnctic ground state. Without quantum di- 
vergence, it is not necessary to modify the dispersion 
relations (2.26a) into the temperature-dependent form 
(2.31). While the thermodynamics should be modified, 
the quantum mechanics may be left as it is. 

Such an idea leads to the Bogoliubov transformation 
free from temperature replacing Eq. (2.37) by 71 (A;) = 0, 
that is, 



tanh2^fe 



2y/Ss cos ak 
S + s 



(2.39) 



The ground-state energy and the dispersion relations are 
simply given by 

Eg = -2NJSs + Ei; LO^=Lot{k), (2.40) 

within the up-to-0(6'^) treatment and by 

Eg = -2NJSs + Ei+Eo; uj^ = ujf{k) + uj^{k), 

(2.41) 

in the up-to-0(S'") treatment. They are nothing but the 
T = findings in the Takahashi scheme. 

At finite temperatures we replace Qj,Qfe and /3^/3fe in 
the spin-wave Hamiltonian (2.23) by 



n=F — 



(2.42) 



n~ ,n+=0 



where Pfe(n^,n+) is the probability of nr ferromagnetic 
and n+ antiferromagnetic spin waves appearing in the 
fc-momentum state and satisfies 

E Pk{n-,n+) = l, (2.43) 

n~ ,n+ 

for all A:'s. Then the free energy is written as 

k n~ ,n+ T=± 

+kBTj2 E Pk{n-,n+)lnPkin-,n+). (2.44) 



We minimize the free energy with respect to Pfe(n ,n+) 
enforcing a condition 



(S^-<)t + 25=(:5^-< :)t 
= S + s- 



N 



k T=± 



i^k 



(2.45) 



as well as the trivial constraints (2.43). In the second-side 
compact expression, the normal ordering is taken with re- 
spect to both operators a and /3. Equation (2.45) claims 
that the thermal fluctuation {S + •s)X]fe(^fc ~^ ''^k)/^k 
should cancel the full, or classical, Neel order (S + s)N 
rather than the quantum mechanically reduced one {S + 
s — 2S)N. Without consideration of the quantum fluc- 
tuation 26, which is absent from ferromagnets but pe- 
culiar to ferrimagnets, the present scheme breaks even 
the conventional spin-wave achievement at low temper- 
atures. Numerically solving the thermodynamic Bethe- 
Ansatz equations. Takahashi and Yamada [57] suggested 
that the conventional spin-wave theory correctly gives 
the low-temperature leading term of the specific heat. 
Both the Takahashi scheme with Eq. (2.28) and the new 
scheme with Eq. (2.45) indeed keep unchanged the con- 
ventional spin-wave findings 
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c 



^-^C(|) ,i/ 



SS 



t'l' (T^O), 



(2.46) 



where t = k^T j J within the up-to-0(S'^) treatment, 
while t = kBT/-fJ with 7=1 + T/VSs - {S + s)A/Ss in 
the up-to-0(5°) treatment. The conventional spin-wave 
approach gives no quantitative information on the mag- 
netic susceptibility, whereas the modified theory reveals 



xJ 



Ss{S-sf 



(T^O). 



In terms of the optimum distribution functions 



(2.47) 



(2.48) 



g[Jw±-i/(S+s)/wfc]/feBT _ ^ ' 

the free energy at the thermal equilibrium is written as 
F = Eg + u{S + s)N -kBTj2Yl (1 + ^fc) ' (2-49) 

k r=± 

where v is the Lagrange multiplier due to the constraint 
(2.45). 

III. RESULTS 

First we calculate the ground-state energy Eg and the 
antiferromagnetic excitation gap and compare them 
with numerical findings in Table I. At T = 0, the Taka- 
hashi scheme and the new scheme lead to the same re- 
sults both giving 1/ = 0. We are fully convinced that 
the spin- wave treatment better works for larger spins. 
We further learn that the spin-wave approach is bet- 
ter justified with increasing 5/s as well as Ss, which is 
because the quantity S — s fills the role of suppressing 
the divergence in Eq. (2.38). On the other hand, the 
Schwinger-boson approach constantly gives highly pre- 
cise estimates of the low-energy properties. Figure 1 fur- 
ther demonstrates that the Schwinger-boson mean-field 
theory is highly successful in describing the low-lying ex- 
citations. Both the bosonic languages well interpret the 
ferromagnetic excitations, whereas the linear spin waves 



considerably underestimate the antiferromagnetic excita- 
tion energies. The quantum correlation has much effect 
on the antiferromagnetic excitation mode and such an ef- 
fect is well included into the Schwinger-boson calculation 
even at the mean-field level. 

Next we calculate the thermodynamic properties. Fig- 
ure 2 shows the temperature dependence of the specific 
heat. The Schwinger-boson mean-field theory is still 
highly successful at low temperatures, while with increas- 
ing temperature, it rapidly breaks down failing to repro- 
duce the Schottky-type peak. The mean-field order pa- 
rameter il monotonically decreases with increasing tem- 
perature and reaches zero at 



kBT 
J 



S- 



1 



ln(l + f/5) + ln(f + l/s) 



(3.1) 



Above this temperature, Q sticks at zero suggesting no 
antiferromagnetic correlation in the system. The onset 

of the paramagnetic phase at a finite temperature is a 
mean-field artifact and the particular temperature (3.1) 
is an increasing function of 5 and s. The modified spin- 
wave theory based on the Takahashi scheme also fails to 
describe the Schottky peak. Because of the Lagrange 
multiplier v, which turns out a monotonically increasing 
function of temperature, the dispersion relations (2.31) 
lead to endlessly increasing energy and thus nonvanish- 
ing specific heat at high temperatures. Only the modified 
spin-wave theory based on the new scheme succeeds in 
interpreting the Schottky peak. Since the antiferromag- 
netic excitation gap is significantly improved by the in- 
clusion of the 0(5°) correlation, the interacting modified 
spin waves reproduce the location of the Schottky peak 
fairly well. Mixed-spin trimeric chain ferrimagnets have 
recently been synthesized [5] and their low-temperature 
thermal properties were well elucidated by the modified 
spin- wave theory [34]. However, it was unfortunate that 
the additional constraint was imposed on the uniform 
magnetization and therefore the higer-tcmpcrature prop- 
erties were much less illuminated. Controlling the stag- 
gered magnetization instead based on the new scheme, we 
can fully investigate such polymeric chain compounds as 
well. 
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FIG. 1. The Schwinger-boson (SB), linear- modificd-spin- wave (LMSW), perturbational interacting- modified-spin- wave 
(PIMSW), and quantum Monte Carlo (QMC) calculations of the dispersion relations of the ferromagnetic (w^) and anti- 
ferromagnetic (w;^) elementary excitations for the spin-(S, s) ferrimagnetic Heisenberg chains at zero temperature. 
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FIG. 2. The Schwinger-boson (SB), linear-modified-spin-wave (LMSW), perturbational interacting- modified-spin- wave 
(PIMSW), and quantum Monte Carlo (QMC) calculations of the specific heat C as a function of temperature for the spin-(S, s) 
ferrimagnetic Heisenberg chains. The modified spin waves are constructed in two different ways, the Takahashi scheme (Taka- 
hashi) and the new scheme (Yamamoto). 
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FIG. 3. The ferromagnetic (oj^T^q) and antiferromagnetic (uj'j^^g) excitation gaps as functions of temperature for the spin-(S', s) 
ferrimagnetic Heisenberg chains calculated by the Schwinger bosons (SB) and the perturbationally interacting modified spin 
waves (PIMSW) based on the Takahashi scheme. 



In the Schwinger representation and the modified spin- 
wave treatment based on the Takahashi scheme, the en- 
ergy spectrum depends on temperature. Since the low- 
energy band structure is well reflected in the thermal 
behavior and can directly be observed through inelastic- 
neutron-scattering measurements, we investigate the fer- 
romagnetic (i^jT^q) and antiferromagnetic (w^^q) excita- 
tion gaps as functions of temperature in Fig. 3. The 
Schwinger-boson mean-field theory claims that the an- 
tiferromagnetic gap should first decrease and then in- 
crease with increasing temperature, while the modified 
spin-wave theory predicts that the excitation energies 
of both modes should be monotonically increasing func- 
tions of temperature. We find a similar contrast between 
the two languages applied to ladder ferrimagnets [35,58]. 
In the case of Haldane-gap antiferromagnets, both the 
Schwinger-boson and modified-spin- wave [41] findings, 
together with the nonlincar-cr-model calculations [59,60], 
commonly suggest that the Haldane gap is a simply acti- 
vated function of temperature. Extensive measurements 
on spin-1 antiferromagnetic Heisenberg chain compounds 
[61-63] also report that the Haldane massive mode is 
shifted upward with increasing temperature. We en- 
courage neutron-scattering experiments on ferrimagnetic 
chain compounds to solve the present disagreement be- 
tween the Schwinger-boson and modified-spin-wave cal- 
culations of the antiferromagnetic excitation gap as a 



function of temperature. 

Figure 4 shows the temperature dependence of the 
magnetic susceptibility-temperature product, which elu- 
cidates ferromagnetic and antiferromagnetic features co- 
existing in ferrimagnets [25] . xT diverges at low tempera- 
tures in a ferromagnetic fashion but approaches the high- 
temperature paramagnetic behavior showing an antifer- 
romagnetic increase. The modified spin waves much bet- 
ter describe the magnetic behavior than the Schwinger 
bosons. The spin waves modified along with the Taka- 
hashi scheme better work at high temperatures, while 
those along with the new scheme precisely reproduce 
the low-temperature behavior. Both calculations con- 
verge into the paramagnetic behavior xkuT / N {g ^ib)^ = 
[S'(S'+ 1) -I- s{s + l)]/3 at high temperatures, whereas the 
Schwinger-boson mean-field theory again breaks down at 
the particular temperature (3.1). Considering that nu- 
merical tools less work at low temperatures, we realize 
the superiority of the new-scheme-based modified spin- 
wave theory all the more. 

Finally we calculate another type of ferrimagnet in 
order to demonstrate the constant applicability of the 
present new scheme. Figure 5 shows the thermody- 
namic properties of the fcrromagnetic-ferromagnetic- 
antiferromagnetic-antiferromagnetic bond-tetrameric 
spin-i Heisenberg chain. 
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FIG. 4. The Schwinger-boson (SB), linear-modified-spin-wave (LMSW), perturbational interacting-modified-spin-wave 
(PIMSW), and quantum Monte Carlo (QMC) calculations of the susceptibility-temperature product xT as a function of 
temperature for the spin-(S', s) ferrimagnetic Heisenberg chains. The modified spin waves are constructed in two different ways, 
the Takahashi scheme (Takahashi) and the new scheme (Yamamoto). 
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FIG. 5. The linear-modified-spin-wave (LMSW), pertur- 
bational interacting-modified-spin-wave (PIMSW), and quan- 
tum Monte Carlo (QMC) calculations of the specific heat C 
and the susceptibility-temperature product xT as functions 
of temperature for the spin-i bond-tetrameric ferrimagnetic 
Heisenberg chain of Jf = Jaf- The modified spin waves are 
constructed on the new scheme. 

N 

4n-3 ■ + S ' S in-i) 



- Jf(<S'4„_i • 5^4 

n ^ An ' '^4n+l) ; 



(3.2) 



where we have set all the g-factors equal for simplicity. 
The new raodified spin-wave scheme again successfully 
reproduces the Schottky peak of the Specific heat. The 
interacting modified spin waves further interpret the low- 
temperature Shoulder-like structure. The characteris- 
tic minimum of the susceptibility-temperature product is 
unfortunately less reproduced but the calculation again 
correctly gives the paramagnetic susceptibility at suffi- 
ciently high temperatures. A recent experiment [64] on 
a single-crystal sample of Cu(C5H4NCl)2(N3)2 [7], which 
may be described by the Hamiltonian (3.2), has reported 
that the specific heat exhibits a double-peaked structure 
as a function of temperature. There is indeed a possi- 
bility of an additional peak appearing at low tempera- 
tures as the ratio Jf/Jaf moves away from unity [65]. 
However, no parameter assignment has yet succeeded in 
interpreting all the observations consistently. There are 
further chemical attempts to synthesize novel ferrimag- 
nets. Organic ferrimagnets [6,66] are free from magnetic 
anisotropy and thus suitable for analyzing in terms of the 
modified spin waves. 




o QMC (b) 
— SB 

LMSW(Takahashi) 
PIMSW(Takahashi) 
FDIMSW(Takahashi) 



1.0 kjiT/J 2.0 3.0 



FIG. 6. The 

linear-modified-spin-wave (LMSW), perturbational interact- 
ing-modified-spin-wave (PIMSW), fuU-diagonalization inter- 
acting-modified-spin-wave (FDIMSW), and quantum Monte 
Carlo (QMC) calculations of the specific heat C and the sus- 
ceptibility X as functions of temperature for the spin-1 anti- 
ferromagnetic Heisenberg chain. The modified spin waves are 
constructed on the Takahashi scheme. 



IV. SUMMARY AND DISCUSSION 

We have demonstrated the Schwinger-boson mean- 
field representation and the modified spin-wave treat- 
ment of one-dimensional Heisenberg ferrimagnets. The 
Schwinger bosons form an excellent language at low tem- 
peratures but rapidly lose their validity with increasing 
temperature. The modified spin-wave theory is more re- 
liable in totality provided the number of bosons is con- 
trolled without modifying the native energy structure. 
On the other hand, the Schwinger-boson representation 
can be extended to anisotropic systems [67] more rea- 
sonably because it is rotationally invariant in contrast 
to the modified spin-wave theory. While the tempera- 
ture dependence of the anitiferromagnetic excitation gap 
WjLq is left to solve experimentally, we are now convinced 
that the bosonic languages remain effective in low di- 
mensions and may be applied to extensive ferrimagnets 
[68] . Besides ground-state properties and thermodynam- 
ics, quantum spin dynamics [69,70] can be investigated 
through the modified spin- wave scheme. 

We further mention our findings in the antiferromag- 
netic limit with the view of realizing the close relation 
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between the two bosonic languages. We equalize s with 
S and set 2N, the number of spins, equal to L for the 
Hamiltonian (2.1). At S — s, the ground-state sublat- 
tice magnetization (2.38) diverges and therefore the new 
modified spin-wave scheme is no more applicable. We 
have to settle the quantum, as well as thermal, divergence 
inevitably employing the Takahashi scheme. Besides the 
perturbational treatment of Ho, we may consider the full 
diagonalization of Tlx + Tio , where the ground-state en- 
ergy and the dispersion relations are still given by Eqs. 
(2.32) and (2.33), respectively, but with 6i~ satisfying 

7i(^)+7o(fc) =0. (4.1) 

Such an idea applied to ferrimagnets ends in gapped 
ferromagnetic excitations and misreads the low-energy 
physics. The perturbational series-expansion approach 
is highly successful in the case of ferrimagnets [32,33]. 
Focussing our interest on Haldane-gap antiferromagnets, 
we list the bosonic calculations of the ground-state prop- 
erties in Table II. The bosonic languages interpret 
the ground-state correlation very well but underesti- 
mate the Haldane gap considerably. Indeed they can- 
not detect the topological terms responsible for vanish- 
ing gap [72], but they are still qualitatively consistent 
with the nonlinear-(T-model quantum field theory, yield- 
ing the low-temperature limiting behavior cj^^p — Aq cx 
g-Ao/T [41 50] and the large-spin asymptotic behavior 
Aq oc e""'^ [44,46,72]. The Schwinger-boson mean- 
field theory and the fuU-diagonalization interacting mod- 
ified spin-wave treatment give the same estimate of the 
Haldane gap. The Schwinger-boson dispersion relation 
(2.11) indeed coincides analytically with that of the full- 
diagonalization interacting modified spin waves at zero 
temperature. This is interesting but not so surpris- 
ing, because the Holstein-Primakoff bosons (2.19) are 
obtained by replacing both a„-|- {bni) and aj^^ (^Ij.) 

by Y^25 — al^^ttni {'\J'^s — h\^.^hn-[ ) in the transformation 
(2.2). 

Figure 6 shows the thermodynamic calculations for the 
spin-1 antiferromagnetic Heisenberg chain. We learn that 
the Schwinger-boson mean-field theory does not work at 
all for spin-gapped antiferromagnets at finite tempera- 
tures, which is in contrast with its fairly good represen- 
tation of the low-temperature thermodynamics for fer- 
rimagnetic chains. On the other hand, the modified 
spin-wave treatment maintains its validity to a certain 
extent. Indeed the Takahashi scheme still fails to re- 
produce the antiferromagnetic Schottky-type peak of the 
specific heat, but it describes the susceptibility very well 
except for the low-temperature findings attributable to 
the underestimate of the Haldane gap. We may expect 
the modified spin waves to efiiciently depict the dynamic, 
as well as static, susceptibility for extensive spin-gapped 
antiferromagnets including spin ladders [73]. As for the 
thermal properties of one-dimensional antiferromagnets, 
whether spin gapped or not, there is a possibility of a 



fermionic language [74.75], which is in principle compact, 
being superior to any bosonic representation. 

In the case of ferromagnets, the Holstein-Primakoff 
bosons are already diagonal in the momentum space 
[37,56], suggesting no quantiuu fluctuation in the ground 
state, and therefore the present new scheme turns out 
equivalent to the Takahashi scheme. The new-scheme- 
based modified spin-wave theory is the very method for 
low-dimensional ferrimagnets and is ready for extensive 
explorations. 
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TABLE I. The Schwinger-boson (SB), linear- modified-spin- wave (LMSW), perturbational interacting-modified-spin-wave 
(PIMSW), and numerical diagonalization (Exact) calculations of the ground-state energy Eg and the zero-temperature antifer- 
romagnetic excitation gap Ao for the spin-(S', s) ferrimagnetic Heisenberg chains. 



Approach 


{S,s) 


= (1,|) 


{S,s) 


= (- -) 

V2' 2/ 


iS,s) 


= (|,1) 


EJNJ 


Ao/J 


Eg/NJ 


Ao/J 


Eg/NJ 


Ao/J 


SB 


-1.45525 


1.77804 


-1.96755 


2.84973 


-3.86270 


1.62152 


LMSW 


-1.43646 


1 


-1.95804 


2 


-3.82807 


1 


PIMSW 


-1.46084 


1.67556 


-1.96983 


2.80253 


-3.86758 


1.52139 


Exact 


-1.4541(1) 


1.759(1) 


-1.9672(1) 


2.842(1) 


-3.861(1) 


1.615(5) 



TABLE II. The Schwinger-boson (SB), linear-modified-spin-wave (LMSW), perturbational interacting-modified-spin-wave 

(PIMSW), full-diagonalization interacting-modified-spin-wave (FDIMSW), and quantum Monte Carlo (QMC) [71] calculations 
of the ground-state energy Eg and the lowest excitation gap Ao for the spin-S' antiferromagnetic Heisenberg chains. 



Approach 



5=1 



S = 2 



S = 3 



EJLJ 



Ao/J 



EJLJ 



Ao/J 



EJLJ 



Ao/J 



SB 

LMSW 

PIMSW 

FDIMSW 

QMC 



-1.396148 
-1.361879 
-1.394853 
-1.394617 
-1.401481(4) 



0.08507 
0.07200 
0.07853 
0.08507 
0.41048(6) 



-4.759769 
-4.726749 
-4.759760 
-4.759759 
-4.761249(6) 



0.00684 
0.00626 
0.00655 
0.00684 
0.08917(4) 



-10.1231 
-10.0901 
-10.1231 
-10.1231 
-10.1239(1) 



0.00295 
0.00279 
0.00287 
0.00295 
0.01002(3) 
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